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DESCRIPTION 
BACKGROUND OF THE INVENTION 

Field of the Invention 
The present invention relates generally to a computerized method and 
system for providing a solution to any general (non-convex) minimum-maximum 
problem in constrained optimization in any specific application. In an exemplary 
application, the technique is used to demonstrate optimal design of three- 
dimensional parts for manufacturing, constrained by prescribed tolerances. 

Description of the Related Art 
The present invention is introduced by first setting forth the following 
known construct: 

Given a function y = f(x,c) = C1X1 + c 2 x 2 , where x is a set of independent 
variables x = {xi,x 2 }, Xi and x 2 are subsets of x, c = {ci,c 2 } is a set of functional 
parameters, partitioned into two subsets Ci and c 2 , and y is a dependent variable, it 
is desired to minimize (over x 2 ) the maximum (over Xi) of y, subject to a linear 
constraint Ai 2 Xi+A 2 iX 2 ^ bi 2 , where A !2 , A 2 i are sub-matrices and bi 2 is a vector. 
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This means finding appropriate values for vectors Ci and c 2 , so as to solve: 
Min {c 2 x 2 + Max C1X1}, subject to A )2 Xi +A 2 jx 2 ^ b !2 . 

This problem is, therefore, associated with a linear constraint set and a 
piece-wise linear objective function. Figure 1 graphically depicts the problem, 
wherein are shown axes X\ and x 2 , representing two sets of control/decision 
variables. Variable xi might represent, for example, the tolerances for various 
components to be manufactured, in which a simple 3-D component is associated 
with three tolerance variables, one for each dimension. Variable x 2 then describes 
the costs of manufacturing the components as a function of the components' 
tolerance specifications. 

In this setting, the tighter the values of a component's tolerance variables, 
the better is the degree of precision achieved in the component. But, at the same 
time, it is more costly to manufacture the component. 

The vector C = (ci, c 2 ) captures the sensitivity of the two sets of decisions 
Xi and x 2 that can be used to calculate the overall manufacturing design efficiency. 
The polyhedron P, label 100, describes the set of all feasible (xi, x 2 ) values that 
satisfy all engineering and other design requirements (as specified by the tolerance 
parameters matrix, A = [A !2 , A 2i ]) to manufacture the components. 

The general problem is NP-hard and difficult to resolve in a reasonable 

time. The usual techniques used in attempting to arrive at a global optimum are 

Simulated Annealing, Genetic Algorithm or other Monte Carlo type approaches. 
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However, all of these techniques are slow, cumbersome, and do not 
guarantee a global solution. 

The generalized min-max problem is characterized as a model for 
scenarios in which there are two conflicting objectives and a concave objective 
function providing an effectiveness/efficiency metric. Presently, no computerized 
tool is available that guarantees a global solution to the general min-max problem, 
let alone a computer tool that can find such solution in a reasonable time. 

SUMMARY OF THE INVENTION 

In view of the foregoing, and other, exemplary problems, drawbacks, and 
disadvantages of the conventional methods, it is an exemplary feature of the 
present invention to provide a computer structure and method for providing a 
global optimum for applications of the generalized min-max problem. 

It is another exemplary feature of the present invention to provide a 
method in which this global optimum can be located in a reasonable amount of 
time. 

It is another exemplary feature of the present invention to demonstrate a 
computer tool, as exemplarily implemented specifically in the application for the 
manufacturing of three-dimensional (3D) parts. 
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It is another exemplary feature of the present invention to provide a 
computer tool available to any number of environments that can be described as a 
min-max problem. 

To achieve the above exemplary features and others, in a first exemplary 
aspect of the present invention, described herein is a method for a computerized 
method for providing an optimization solution, including, for a process, wherein 
is defined a linear functional form y = f (X,c), where X comprises a set of 
independent variables X = {xi,...x n }, c comprises a set of functional parameters c 
= {ci,...c n } 5 and y comprises a dependent variable, where the independent 
variables set X is partitioned into two subsets, Xi and X 2 , receiving data for said 
process, minimizing y with respect to Xi, and maximizing y with respect to X 2 , 
subject to a set of constraints, wherein the maximizing y comprises a global 
optimum for the process. 

In a second exemplary aspect of the present invention, described herein is 
an apparatus for calculating a global optimization to a minimum-maximum 
problem, including a first calculator to provide a plurality of minimum values, and 
a second calculator to locate a global optimum value, given the plurality of 
minimum values. 

In a third exemplary aspect of the present invention, also described herein 
is a system including a memory containing data appropriate to a minimum- 
maximum problem and an apparatus further including a first calculator to provide 
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a plurality of minimum values and a second calculator to locate a global optimum 
value, given the plurality of minimum values. 

In a fourth exemplary aspect of the present invention, also described herein 
is a signal-bearing medium tangibly embodying a program of machine-readable 
instructions executable by a digital processing apparatus to perform the above 
method for providing an optimization solution. 

In a fifth exemplary aspect of the present invention, also described herein 
is a business method, including at least one of: executing a computerized 
calculation to find a global maximum in accordance with the above-described 
method; providing a data for the process, the data to be used in the computerized 
calculation for the global optimum; receiving a result from the computerized 
calculation; providing one or more software modules for the computerized 
calculation; and developing one or more software modules for the computerized 
calculation. 

In a sixth exemplary aspect of the present invention, also described herein 
is a computer tool for providing a global solution to a minimum-maximum 
problem, including a linear programming solver to calculate a periphery of a 
polyhedron representing a region of all points that satisfy a linear constraint in a 
minimum-maximum problem. 
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The present invention, therefore, provides a method that determines a 
global solution for a generalized min-max problem and a computer tool that will 
execute the method in a reasonably short time. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other exemplary purposes, aspects and advantages will 
be better understood from the following detailed description of an exemplary 
embodiment of the invention with reference to the drawings, in which: 

Figure 1 shows a feasible region (polyhedron) 100 for a general min-max 
problem; 

Figure 2 shows an overview 200 of an exemplary system in which the 
present invention is applied to manufactured parts tolerance design; 

Figure 3 shows an exemplary flow chart 300 for the min-max optimization 
step 204 of Figure 2; 

Figure 4 shows an exemplary flow chart 400 for obtaining the min-max 
solution step 303 of Figure 3; 

Figure 5 shows an exemplary flow chart 500 for the min-max problem 
optimization step 401 of Figure 4; 

Figure 6 shows that the objective function to be minimized in the 
resolution of the min-max problem is a piece-wise linear, concave function, g(x 2 ) 
600; 

YOR920030256US1 

6 



Figure 7 shows the set of allowable solutions, D = (x 2 , g(x 2 )), and the 
polyhedron, conv(D), the convex hull of the points in D, as the feasible region 
associated with the reformulated problem; 

Figure 8 exemplarily shows a 1 -polar cut for conv(D); 

Figure 9 illustrates an exemplary hardware/information handling system 
900 for incorporating the present invention therein; 

Figure 10 illustrates a signal bearing medium 1000 (e.g., storage medium) 
for storing steps of a program of a method according to the present invention; and 

Figure 1 1 shows an exemplary block diagram of the software modules for 
a computer tool that executes and embodies the method described herein. 



DETAILED DESCRIPTION OF AN EXEMPLARY 
EMBODIMENT OF THE INVENTION 

Referring now to the drawings, and more particularly to Figures 1-1 1 , an 

exemplary embodiment of the computer tool of the present invention will now be 

1 5 discussed in the min-max application of maximizing the tolerance of 

manufactured parts while minimizing the cost for these parts. One of ordinary 

skill in the art, taking the following discussion as a whole, will readily recognize 

that the technique used in this computer tool can easily be adapted to other 

applications, such as product distribution, supplier/vendor availability and 

20 distribution, optimal procurement activities, securities portfolio management, etc. 
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The tolerancing of 3-D parts in manufacturing design allows the parts 
designer to directly associate the tolerance data with the native 3-D geometry. 
Former practices, based on usage of drawings, require the generation of these 
drawings and lead to lack of integration between the technological data 
(tolerances) and geometrical 3-D data. 

One challenge in coming up with the tolerancing scheme associated with 
any part is that it must be consistent with the geometry for this part and that it 
must follow the rules set by the international tolerancing standards (ANSI, ISO, 
etc.). In the design of parts, there are two conflicting goals: 

1) The tolerance values must not be too loose, in order to guarantee that 
the specifications for the assembly will be reached; and 

2) The tolerance values must not be too tight, in order to reduce the 
manufacturing cost of each part. 

In such a problem setting, the computer tool of the present invention aids a 
parts designer in determining the optimum tolerancing scheme that minimizes the 
manufacturing costs for all the parts by accommodating for maximum allowable 
tolerances (e.g., to meet the assembly specifications). 

Underlying the computer tool of the present invention is a novel technique 

that yields a global optimum in the general case of non-convex problem for the set 

of problems described above. The present invention is cognizant of the 

aforementioned functional construct shown in Figure 1 . Moreover, the present 
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invention builds upon this known functional construct, but references this known 
construct to impose upon it novel problems, constraints, and desiderata, of the 
following illustrative type. 

A typical conventional approach to solving this nonlinear problem would 
be that of parameterizing x 2 and evaluating the piecewise linear function g(x) at 
every choice of x 2 . The present invention also uses an iterative procedure to solve 
the min-max problem, but it implements an iterative process in which local 
minima are first found and then the maximal member among the local minima is 
located, thereby providing a global optimum. 

Figure 2 shows an exemplary overview 200 of the information flow in the 
computer tool of the present invention, as applied to be a process for determining 
optimal parts design. 

In Figure 2, database 201 contains tolerance data for the parts, and 
database 202 contains data for costs of manufacturing these parts. In step 203, the 
data for tolerance and costs are provided as inputs into the present invention, and 
in step 204, the tolerance design problem is solved using a min-max optimization 
model discussed below to provide an output 205 being the optimal parts design. 

Figure 3 shows an exemplary process 300 to implement the min-max 

optimization model (e.g., reference step 204 of Figure 2). Step 301 is the entry 

into the process, so that, in step 302, the data from the two databases 201,202 are 

fitted into data structures to populate the min-max model. That is, the data in 
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these two databases 201, 202 would be converted into appropriate data structures, 
one data structure for cost parameters (c) and one for tolerance parameters (A). 

The costs parameters, c, define the sensitivity of the two sets of tolerance 
decisions Xi and x 2 that can be used to calculate the overall efficiency of 
manufacturing. The tolerance parameters matrix A specifies the engineering and 
other design requirements that guide the process of obtaining optimal tolerance 
specifications for all the components. 

In step 303, the min-max problem is solved to provide the data to generate 
the design reports in step 304. It is noted that a number of algorithms could be 
implemented in software format for this step 303, including a number of known 
algorithms. 

However, as initially mentioned above, the problem with using these 
known algorithms is that they provide, at best, a local optimum (e.g., they do not 
guarantee a global optimum). In contrast, the method of the present invention, 
shown in overview in Figures 4 and 5, uses an iterative optimization looping 
technique that further provides a global optimum. 

Thus, Figure 4 shows an exemplary process 400 for step 303 (e.g., solving 
the min-max problem instance) of Figure 3. However, as an overview of the 
novel process of the present invention, it would be helpful to return to Figure 1 to 
graphically explain the mathematical process. 
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As mentioned above, the polyhedron P in Figure 1 describes the set of all 
feasible (xi, x 2 ) values that satisfy all engineering and other design requirements 
(as specified by the tolerance parameters matrix, A = [Ai 2 , A 2 i]) to manufacture 
the components. Therefore, the min-max solution lies somewhere in the 
polyhedron P. The technique of the present invention recognizes that the solution 
actually lies on the periphery of polyhedron P and that the global solution to the 
min-max problem can be found by "sliding" sensitivity vector C across the 
polyhedron P to determine which corner point of the polyhedron provides the 
min-max global solution. 

Thus, important contributions of the present invention include a 
recognition that a computer tool would be able to "construct" the polyhedron P 
from the data in the tolerance and cost databases. The tool can then systematically 
test the periphery of the polyhedron P, and more specifically, the corners of the 
periphery, by constructing "1 -polar cuts" by interconnecting corner points of the 
polyhedron P and using sensitivity vector C to zoom in on the region of the 
1 -polar cut the solution lies in. 

In view of this graphical description and returning to Figure 4, in step 

303.1, the objective function for total design cost is calculated. This step specifies 

the objective function in terms of a linear equation defined in terms of the cost 

vector c 5 and the tolerance decisions x, expressed as (c 2 x 2 + CiXi). In step 303.2, 

the relationship between design cost and tolerances is introduced, and in step 
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303.3, any other structural interrelationships between design parts are introduced. 
In these two steps, the engineering and other design requirements, specified by 
data items in A, are enforced. 

Finally, in step 401, the min-max problem is then optimized (e.g., a global 
optimum is found), using the exemplary looping technique 500 shown in Figure 5. 

In step 401.1 of Figure 5, the design costs are minimized by solving the LP 
relaxation. An LP relaxation of the min-max problem (potentially a non-linear 
mathematical program) refers to the linear programming representation that 
results by dropping the non-linearities present in the model formulation. There are 
standard commercially available linear programming (LP) solvers such as 
CPLEX, OSL, XPRESS, that can be used to solve such an LP relaxation. 

In step 401 .2, it is determined whether the solution of step 401 . 1 is 
optimal, and, if so, the process is stopped in step 401 .4. If the solution is not 
optimal, then in step 401.3, a 1 -polar cut that separates the incumbent solution 
from the min-max polytope is found, and the process returns to step 401 . 1 . This 
looping continues until the global optimum is located. 

A "1 -polar cut" is a mathematical construct based on polarity theory (as 
described in the text book titled Integer Programming and Combinatorial 
Optimization by Nemhauser and Wolsey (John Wiley and Sons, New York, 1988) 
that can be computationally generated to obtain an additional chracterization for 
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conv(D) as illustrated in Figure 8, so that local optima can be discarded by this 
"cut" to proceed towards a global optimum. 

Figures 6 through 8 illustrate in greater detail the mathematical foundation 
underlying the computerized method implemented by the present invention. 

More specifically, for natural numbers k = 1 Pmaster(u k ) is defined as: 
Pmaster(uk) = Min{c 2 x 2 + u k (bi2-A2iX 2 ): Ai 2 Xi + A 2 ix 2 ^ bi 2 }. 

Psub(x 2k ) is defined as the LP problem: 

Psub(x 2k ) = Max{ciX! + c 2 x 2k : A !2 xi ^ bi 2 - A 2 ix 2k }, with associated 
optimal dual multipliers u k . 

Let zM k be the optimal objective function value of program Pmaster(u k ) 
and zS k be the same for program Psub(x 2k ) for k = 0, 1 

With these definitions, procedure Local, discussed below, is used to begin 
to find an optimal solution to the min-max problem. Procedure Local assumes the 
existence of a known extreme point solution in the polyhedron P. This is not a 
strong assumption, as a candidate starting solution can always be generated by 
solving the linear program Pmaster(O). 
Procedure Local 

Given an initial feasible extreme point solution g(x) = (g(xi), g(x 2 )) in P, 

or, equivalent symbology, x g = (xi g , x 2g ) in P, the steps of the procedure Local are: 

Step 0 : Set k = 0, z M 0 = Cix ]g + c 2 x 2g , and x 20 = x 2g . 

Step 1 : Solve Psub(x 2k ) to obtain u k and z s k . 
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Step 2 : Solve Pmaster(u k ) to obtain x 2k+ i and z M k+ i. 

If z M k+ i ^ z M k then STOP; else, set k to k+1, goto Step 1 . 

Steps 1 and 2 above can both be executed using an LP solver. In terms of 
the exemplary application for parts manufacturing, in step 1, at given tolerance 
costs, x 2k , the k-th step results in an intermediate solution associated with the 
design efficiency value z s k , and a design efficiency improvement direction, u k , as a 
recommendation to proceed towards the local optimum. In step 2 above, in the 
suggested improvement direction, u k , the k-th step finds a better tolerance 
alternative, x 2k+ i, that is guaranteed to improve upon the incumbent design 
efficiency value. 

The result of Procedure Local is a set of tolerance decisions associated 
with a deisgn efficiency value, z M k , that may be improved further. 

It can be shown that procedure Local converges to an extreme point 
solution of P in finitely many steps. However, the solution to which this 
procedure converges can at best be a local optimum, as can be shown using the 
following result: g(x) is concave for x in P 2 (e.g., see Figure 6). 

Thus, the global optimizer to problem Pmin-max , x 2 *, lies on the relative 
boundary of P 2 . What this says is that, if any gradient-based search is used to 
solve the min-max problem Pmin-max, as in Figure 6, then the solution may be 
trapped in a local optimum. 
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Therefore, parameterizing x 2 and evaluating g(x 2 ) (e.g., equivalent to 
solving an LP) may not be a viable approach to yield the global optimal solution 
to Pmin-max. In fact, this "parameterization approach" may not even be able to 
improve on a given non-optimal initial solution. 

In the present invention, in order to overcome this limitation, the min-max 
problem is reformulated in a manner that permits one to solve the min-max 
problem to global optimality. This new formulation, in accordance with the 
present invention, of the min-max problem solves the following problem: 

min{ciXi + c 2 x 2 : (xi,x 2 ) belong to conv(D)}, and this problem can be 
solved as an LP via successive refinement of the feasible region. The 
reformulated problem is shown graphically in Figure 7. Here, D 700 represents 
the set of all local minima that can be generated by solving an LP relaxation of the 
min-max problem. Conv(D) 701 represents the convex hull of all such local 
minima. 

Following this reformulation, a global optimum solution is formulated 
below (e.g., see Figure 8): 
Procedure Global 

Given A, bi 2 , c, and the definitions of P and D: 

Step 0 : Solve problem Max {c x : A x ^ b !2 } to obtain an optimal 
solution x M , and set v to x M + c. Set k = 0. 
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Step 1 : Solve problem Min {c x : x in P} to obtain an optimal 
solution x m with the associated objective function value Z k . 

Step 2 : If x m belongs to D, then STOP. x m solves Pmin-max. Else, 
solve the polar problem Pi(conv(D), v 5 x m ) to obtain a polar cut {s(x-v) ^ 1 }. Set P 
to P A {x:s(x-v)^l } and k to k+1, and goto Step 1 . 

It can be shown (e.g., referring to polarity theory literature, such as the text 
book titled Integer Programming and Combinatorial Optimization by Nemhauser 
and Wolsey (John Wiley and Sons, New York, 1988, the contents of which are 
incorporated herein by reference) that the polytope P (defined as a bounded 
polyhedron) shrinks "strictly" in Step 2 of each iteration in such a way that in each 
iteration a new optimal solution with a strictly better objective function value is 
generated. This monotonic behavior leads to the finite convergence of the above 
described Procedure Global, as shown in the iterative process shown in Figure 5. 

Therefore, the present invention presents a computerized tool that provides 
a global optimum in any application having the min-max problem characteristic. 
It is exemplarily demonstrated in the scenario of determining the optimum 
tolerancing scheme that minimizes the manufacturing costs for all the parts by 
accommodating for maximum allowable tolerances (e.g., to meet the assembly 
specifications). However, one of ordinary skill in the art would readily recognize, 
after having read the present application and taking the concepts herein as a 
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whole, that the computer tool discussed above would be applicable in many other 
applications. 

As non-limiting examples, the present invention might be adapted for a 
procurement problem in which the goal is to determine the mix of suppliers, given 
data for quantity discounts. In this application, variables xi represent the 
quantities for a given product to be procured from suppliers in different time 
periods. Variables x 2 then describe the supplier costs of procuring as a function of 
the quantities supplied. 

Another scenario might be the optimal portfolio selection, in which 
variables xi represent the quantities for a given security to be held in a portfolio 
over different time periods. Variables x 2 then describe the risks of shortfalls as a 
function of the quantities of securities held. 

Yet another scenario might be a health care problem involving 
survival-time classification of breast cancer patients, in which the goal is to 
minimize the m-minima of k-linear functions. 
Exemplary Hardware Implementation 

Figure 9 illustrates a typical hardware configuration of an information 
handling/computer system in accordance with the invention and which preferably 
has at least one processor or central processing unit (CPU) 91 1. 

The CPUs 91 1 are interconnected via a system bus 912 to a random access 

memory (RAM) 914, read-only memory (ROM) 916, input/output (I/O) adapter 
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918 (for connecting peripheral devices such as disk units 921 and tape drives 940 
to the bus 912), user interface adapter 922 (for connecting a keyboard 924, mouse 
926, speaker 928, microphone 932, and/or other user interface device to the bus 
912), a communication adapter 934 for connecting an information handling 
system to a data processing network, the Internet, an Intranet, a personal area 
network (PAN), etc., and a display adapter 936 for connecting the bus 912 to a 
display device 938 and/or printer 939 (e.g., a digital printer or the like). 

In addition to the hardware/software environment described above, a 
different aspect of the invention includes a computer-implemented method for 
performing the above method. As an example, this method may be implemented 
in the particular environment discussed above. 

Such a method may be implemented, for example, by operating a 
computer, as embodied by a digital data processing apparatus, to execute a 
sequence of machine-readable instructions. These instructions may reside in 
various types of signal-bearing media. 

Thus, this aspect of the present invention is directed to a programmed 
product, comprising signal-bearing media tangibly embodying a program of 
machine-readable instructions executable by a digital data processor incorporating 
the CPU 91 1 and hardware above, to perform the method of the invention. 

This signal-bearing media may include, for example, a RAM contained 

within the CPU 91 1, as represented by the fast-access storage for example. 
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Alternatively, the instructions may be contained in another signal-bearing media, 
such as a magnetic data storage diskette 1000 (Figure 10), directly or indirectly 
accessible by the CPU 911. 

Whether contained in the diskette 1000, the computer/CPU 91 1, or 
elsewhere, the instructions may be stored on a variety of machine-readable data 
storage media, such as DASD storage (e.g., a conventional "hard drive" or a RAID 
array), magnetic tape, electronic read-only memory (e.g., ROM, EPROM, or 
EEPROM), an optical storage device (e.g. CD-ROM, WORM, DVD, digital 
optical tape, etc.), paper "punch" cards, or other suitable signal-bearing media 
including transmission media such as digital and analog and communication links 
and wireless. In an illustrative embodiment of the invention, the 
machine-readable instructions may comprise software object code. 

In yet another aspect of the present invention, Figure 1 1 shows an 
exemplary block diagram 1 100 of the software modules for a computer program 
that might be used in the exemplary system 900. 

In Figure 1 1, the user would interface with the computer program through 

the graphic user interface (GUI) 1 101 for controlling the execution and for 

viewing results and/or otherwise providing the results to be stored in a data file, 

printed out, or exported into another application program. Memory interface 1 102 

allows the database data to be accessed and intermediate and final results to be 

stored in memory. Calculator module 1 103 includes submodules to convert the 
YOR920030256US1 

19 



database data into the appropriate data structure formats, to execute the local 
subroutine, and to execute the global subroutine. 

In yet another aspect of the present invention, one of ordinary skill in the 
art, taking the discussion above as a whole, would readily recognize that this 
computer tool could also be the basis of a business method or service that relies 
upon some aspect of this tool. Thus, for example, a service might be implemented 
in the form of a consulting enterprise having the function of reading a database 
and executing the min-max solution using the computer tool taught in the present 
invention. This service might even be based on using the Internet as the means of 
communication between the enterprise and the client and the client's database. 

Another possible service based on the present invention might be 
implemented in the form of a software enterprise having the function of adapting 
the computer tool for specific applications, such as those mentioned above. 

As will be readily recognized by one of ordinary skill in the art, this 
invention is well suited for modeling and resolving linear optimization problems 
where there are two conflicting objectives, and there is a well-understood notion 
of an effectiveness/efficency metric (to be used as an objective function) that is 
concave. 

While the invention has been described in terms of exemplary 

embodiments, those skilled in the art will recognize that the invention can be 

practiced with modification within the spirit and scope of the appended claims. 
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Further, it is noted that Applicants' intent is to encompass equivalents of 
all claim elements, even if amended later during prosecution. 
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